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Abstract 

We study integrable models for electrons in metals when the single particle spec- 
trum is discrete. The electron-electron interactions are BCS-like pairing, Coulomb 
repulsion, and spin exchange coupling. These couplings are, in general, nonuniform 
in the sense that they depend on the levels occupied by the interacting electrons. By 
using the realization of spin 1/2-operators in terms of electrons the models describe 
spin 1/2 models with nonuniform long range interactions and external magnetic 
field. The integrability and the exact solution arise since the model Hamiltonians 
can be constructed in terms of Gaudin models. Uniform pairing and the resulting 
orthodox model correspond to an isotropic limit of the Gaudin Hamiltonians. We 
discuss possible applications of this model to a single grain and to a system of few 
interacting grains. 



1 Introduction 



Recent advances in nanophysics allowed experimental investigations on iso- 
lated quantum dots and metallic nanoparticlesfl]. Because of their small size, 
such systems have a discrete spectrum, which can be resolved through I- 
V characteristic measurements. Besides potential technological applications, 
spectral analysis provides insight into the nature of the interactions in met- 
als, comparable to the spectroscopy-based knowledge about the interactions 
in atoms and molecules [2]. This motivates the broad recent interest in var- 
ious physical contexts: the spin orbit interaction was studied in mesoscopic 
Al-Au systems[3]; various ferromagnetic properties were analyzed in experi- 
ments on Co grains[4]; the experiments on Al samples opened up debates on 
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the crossover behavior of superconducting fluctuations, when the size of the 
sample is progressively reduced. The peculiarity of small grains and quantum 
dots is their small capacitance, which fixes the number of electrons in the 
system. In contrast, the characterization of many physical phenomena like 
e.g. phase transitions relies on quantities defined for infinite systems in the 
grand-canonical description. Therefore, the fingerprints of these phenomena, 
i.e. their characterizing physical quantities, must be identified at first [5]. Then, 
techniques suitable for systems in the canonical ensemble must be applied. Fi- 
nally, for such small systems, approximations are difficult to control because 
of the presence of strong finite-size fluctuations and then, exact results are 
extremely valuable. 

Many efforts have been devoted to the construction of theoretical models de- 
scribing the physics of small metallic grains and quantum dots [6-9]. It has 
been argued that, due to the underlying disorder, the leading part of the 
electron-electron interaction reduces the Hamiltonian to a simple "universal" 
form (i.e. independent of the geometry of the system and the realization of the 
disorder), containing only the electrostatic interaction E c , a BCS-like coupling 
g, and a spin-spin coupling J 

H v = Y.^io + E C N 2 - JS 2 - gK + K~, (1) 

where the £j form the effective single-particle spectrum, rii a is the electron 
number operator, S and iV are the total spin and number of electrons respec- 
tively. K~ = J2i c ii c ip K + = are the pair annihilation and creation oper- 
ators. We point out that the couplings in (1) are independent of the involved 
single-particle levels; in this sense, they are uniform. Non-universal corrections 
to Hamiltonian (1) are of relative order 5/ETh, Exh being the Thouless energy 
and S the average level spacing. [6,7] 

Hamiltonian (1) has been used to attack a variety of physical phenomena [10]: 
equilibrium transport near degeneracy have been thoroughly investigated [6]; 
recently, Kurland et al. have shown that the paramagnetic phase is favored by 
the presence of disorder in the single-particle level spacing of the £j [7] . 
For g — 0, the Hamiltonian (1) describes a normal metallic grain with spin- 
spin interaction. Its eigenstates are certain linear combinations of Fock states 
(eigenstates of n^), such to be eigenstates of S 2 as well; they are obtained ap- 
plying the angular momentum lowering operator to the highest weight state. 
For g 7^ the exact solution of H v was found by Richardson and Sherman 
(RS) [11] using techniques close in spirit to the coordinate Bethe Ansatz (BA). 
By expressing Hjj in terms of commuting integrals of motion in Refs. [12,13] 
(see also Ref. [14]), it was demonstrated that this Hamiltonian is also inte- 
grable. In fact, the exact BA eigenstates of Hjj have been successfully employed 
in the recent literature [5,15]. 
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Non-equilibrium transport experiments [16] show that the corrections 0(5/ Etk) 
to Hamiltonian (1) become very important for ultra-small samples. The pres- 
ence of resonance clusters in the tunneling conductivity through Al grains are 
experimental evidence of fluctuations in the electrostatic interaction; in quan- 
tum dots the universal Hamiltonian (or Constant Interaction model) predicts 
a "bimodal" distribution of the tunneling-peak spacing, while a distribution 
corresponding to a Gaussian Orthogonal Ensemble (GOE) is observed [6]. Fur- 
ther, it has been argued that fluctuations in the spin pairing should influence 
the threshold for the mesoscopic Stoner instability, which could be observed 
in future experiments [8]. The presence of fluctuations leads to nonuniform 
couplings in the Hamiltonian, i.e. couplings whose strength depends on the 
two interacting levels. Not all nonuniform couplings yield an exactly solvable 
Hamiltonian; of course, demanding solvability results in restrictions on the 
admitted form of the interactions. 

The aim of the present work is to explore the possible forms of the couplings 
compatible with the exact solvability of the corresponding models. We could 
not find [17] solvable generalizations of Hjj starting from a quadratic bosonic 
model, as done by RS. The strategy we adopt, instead, consists in generalizing 
the procedure of Ref. [12], namely constructing the Hamiltonian of the system 
in terms of Gaudin Hamiltonians. By means of the integrability and the exact 
solvability of the latter the exact solution of the model is obtained. For uniform 
couplings the integrals of motion are tightly related to the isotropic Gaudin 
Hamiltonians; for nonuniform couplings to the anisotropic Gaudin Hamilto- 
nians. The models arising from this procedure include certain nonuniform 
interactions, which could account for corrections to the universal Hamiltonian 
due to e.g. disorder fluctuations. Namely: 

H — ^ e ia n io + Uij UiUj — Qij c\^c\^C^C^ — JijSi • Sj, (2) 
icr ij ij ij 

where U^, and are the nonuniform generalizations of E c , g and J in 
Eq.(l), respectively. This class of models can also be interpreted as describing 
systems of several coupled grains. Then, the coupling constants necessarily 
have to differ significantly when describing intra- and inter-grain interaction 
respectively. Other possible applications may be found in nuclear physics, 
where the nonuniformity of the couplings is notable. 

This article is structured as follows: In section 2, we review how to derive 
the form of the Hamiltonian describing small grains, including the effects of 
disorder. Then we will rewrite the Hamiltonian in terms of spin and charge 
realizations of su(2). This will shed light on its structure, and lead to use- 
ful simplifications in the diagonalization procedure. In section 3, we sketch 
the solution of the universal part of the Hamiltonian and present the isotropic 
Gaudin Hamiltonians and their connection with the integrals of motion. In sec- 
tion 4, we construct integrable Hamiltonians by means of anisotropic Gaudin 
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Hamiltonians. This class of Hamiltonians goes beyond the universal model 
with constant couplings. In general it is quadrilinear in the pseudo-spin oper- 
ators. The tri- and quadrilinear terms can be eliminated by a suitable choice of 
parameters. In section 5, we discuss regimes of the interactions for exemplary 
choices of the parameters. Section 6 is dedicated to concluding remarks. In 
the appendix, we review some exact results obtained by Gaudin. 



2 The model Hamiltonian 



In 2nd quantization, the Hamiltonian for N particles in a confining potential 
V\ = J2 a V(r a ) and a (spin-independent) electron-electron interaction % = 
lEa^bv(r a -r b ), a, be {!,..., N} is 



# = E|| + ^ + ^ = E eiama + E E M%clc) T c lT c k(J , (3) 

a Zm ia ijkl err 

M l k \ = \j dvdv'<t>*{v)<t>*{v')v{v - r%(r')0 fc (r), 

where Si a and 0j(r)xj(cr), i G {1, Q} := / are the eigenvalues and eigenfunc- 
tions of the single particle Hamiltonian (x(cr) is the spinor basis); cj CT , c ia are 
electronic creation and annihilation operators. The disorder average reduces 
the leading part M of the matrix elements M kl to be "uniform-diagonal" ; i.e. 
only elements corresponding to pairwise equal indices are non-zero and con- 
stant: Mil := a, Mj{ := (3 and Mg := 7 for i ^ j, M\\ := uj. We can rewrite 
the interaction matrix as 



M ki = &ik$ji( l - $ij) a + 8u$jk(l - $ij)P + $ijhi(l - kk)l ( 4 ) 

the (non-universal) fluctuations SM k \ are such that SM kl SM kl = 0(S 2 / E^ h ), 
i.e. the corrections are of relative order S/Exh- Inserting expression (4) in the 
Hamiltonian (3), we obtain H = Hu + 5H, where 



icr i^j OT 

+EE cj 



anio-Tijr + (3c icr Cj T c iT Cj a + ^yc i(T c iT Cj T Cj cr 
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+2 (uj - a - (5 - 7) Nj 



pairs 



5h=y: 



'it 



^4 5^71^1 



+ 2 ZX 1 _ MhX 1 - S ik 8ji)(l - 5 u 5 jk )M%cl a c} T c lT c k(7 , 

ijkl °"r 

where rij = X^^mt is the number of electrons in the i-th level-pair, a = —a, 
Npairs = J2i n i] n ii is the number of pairs, and e ia = e ia — a + 2(5 + SM^. 
Introducing the spin operators 

^ = SW > S 7 = { S t)' . ^ = ^(n it -n iA ), (5) 



we have = |rij + 5f, = |rij — S?, and then: 



(a - ^niUj - pSi-Sj + 27cJ T cJ iCji c iT 



(6) 



+2 (u - a - - 7) iVpoir, 



£ 4 5M>» T n, 



I 



(7) 



f/'W err 



From we recover, up to a constant of motion, the universal Hamiltonian 
(1), defining E c = ct — /3/2, J = (3 and <? = —27. This constant of motion is 
relevant to determine the ground state properties. 

The Hamiltonian studied within the present paper is the sum of the uniform 
part (6) and the first line of the nonuniform correction (7): 



a — 2QSI ^2 JijSi ■ Sj 

i ij 



(8) 



where we defined 2£j = — e^, 2^ = £jj + £j|. It can be written in a more 
perspicuous way introducing the operators 

Kf = c] T 4, KJ = (K+) \ K] = \ (n, T + n n - 1) , (9) 
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which are generators of the charge su(2): Kf,K z = ±5{jKf, Kf,K 
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25ijK*. The charge su(2) is orthogonal to the spin su(2), i.e. their genera- 
tors have the property S^K\ = 0. In terms of the spin and charge su{2) 
operators, the Hamiltonian (8) can be finally written as 

H = H S + H K + E (10) 
H s = Y J 2QSt-Y.-h^-S ] , (11) 

i ij 

H K = Y: 2 Vi K* + 4 £ XJijKlK] - l -Y. 9ij{KtK~ + K'Kf), (12) 

i ij ij 

where 2^ = 2^ + 4 Y,j Uij - ga and E = & + E»j *V 

In what follows, we refer to H$ as spin Hamiltonian, and to Hk as charge 

Hamiltonian. They belong to the class of integrable models we find in section 

4. 

We note that [H, S z ] = [H, K z ] = 0. Hence we can diagonalize the Hamiltonian 
in the subspace of fixed number of pairs N K and up spins N s . Because of the 
orthogonality of the spin and charge su(2), it follows that [H^, H$] = and in 
particular [Hs, Kf] = [Hk, Sf] = 0. Thus H $ and Hk are degenerate respect 
to the action of Kf and Sf respectively. Therefore, the Hilbert space H can 
be split into Hs and Hk created by the action of N s and N K Kf on the 
vacua |0) s and |0) x respectively. The spin vacuum |0) 5 = ||, •••,!) is the 
state in which all electrons of the singly occupied level-pairs have spin down; 
the charge vacuum \Q) K is the state in Hk with all level-pairs empty. Namely, 
Hk consists of both the doubly occupied and empty level-pairs which we label 
with the index set Ik {\Ik\ = ^k), and Hs consists of the singly occupied 
levels-pairs, which we label by the index set I s = I \ Ik (I-Tsl = &s)- This 
factorization of the Hilbert space was called the "blocking of singly occupied 
levels" when considering the charge Hamiltonian only [11], and it reduces the 
problem of finding the spectrum of H to the diagonalization of Hs and Hk in 
Hs and Hk, respectively. 



3 Uniform couplings 



In this section, we recast the known results [5, 11, 12] for the universal Hamil- 
tonian in a uniform magnetic field in the present frame. In this case, having 
decomposed the universal model (1) as in formula (10), the spin and charge 
Hamiltonians are 



H S = 2(S Z - JS\ 
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H K = Y J ^Kt- 9 -Y.(KtKj + K~K^)+AE c lr£K-t\ , 

i ij \ i / 



(13) 



and Eq = + QrEc, while the parameters are ( = hbB and 2r/j = 2^ + 
AE c Vt K —g- The eigenstates |^)^ = (S + ) Ns |0) 5 of H s are common eigenstates 
of the total spin S 2 and S z with eigenvalues E Ns = 2fi B B(N s —il s /^) — JS(S+ 
!)• 

H K can be diagonalized through coordinate [11] or algebraic [18] BA. Due to 
its integrability, it suffices to diagonalize its integrals of motion Tj. Indeed, H K 
can be written as 



H K = J2 2r li T i + (4E C + g)J2 T i T i + const., (14) 

i ij 

niK^Kt-gZ?^ + (15) 

j±i '» 'Ij 

where Sj are the isotropic Gaudin Hamiltonians (see also the appendix). With 
a straightforward generalization of Gaudin's results[19,20], one recovers eigen- 
states and eigenvalues of Tj, and from these the eigenvalues E Nk of H K 



a i E a - i^i 

e Nk = J2 e « + e c (2n k - n K f . (16) 

a 

The quasi-energies E a are solutions of the RS equations 

g j ^j - &a p &p - &a 



Here and in the rest of the paper, a primed sum means that coincidences of 
indices are avoided, e.g. J2' a f( a ,f3) '■= J2 ^ f( a ,P)- 



4 nonuniform couplings 



In this section we consider the following class of models, which includes Hamil- 
tonian (10), 



H K = £ 2 Vi Kt -W g, 3 {KlKj + K-K+) + 4 £ U^K] (18) 

i ij ij 
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(19) 



We prove the integrability of this class for a particular choice of the coupling 
constants (Eq. (20) below), and find its exact solution. The proof of inte- 
grability is constructive and performed for the spin and charge Hamiltonian 
separately (see the discussion at the end of section 2). We summarize the 
results of subsections A and B. 

We find that H K and H s are integrable if the couplings are 

pu(r)i - rjj) 



9a 



sinh[p(-Uj — Uj)} 



AUij = A + pu(r)i - rjj) coth[p(w; - Uj)] , 



> for i 7^ j 



smh[q(vi - Vj)} ' 



Jii = ~qv(Ci - 0) coth [q(vi - vA], 



(20) 



where the parameters A, u, v, Ui, r/i and Vi are arbitrary real numbers, while p 
and q are imaginary or real. A quantitative discussion of the couplings for an 
exemplary and physically reasonable choice of parameters, is done in section 5. 
We remark that no analog to A occurs for the spin term. The reason is that we 
can express Hg even without the bilinear term in the integrals of motion, and 
that J2ij r i(S)Tj(S) = (S z ) 2 is an integral of motion due to a symmetry (the 
S^-symmetry of the Hamiltonian) which is not related with integrability. For 
H K) instead, this constant is needed to keep the electron-electron Coulomb 
interaction repulsive. 

For every partition of the set / in distinct subsets Is, Ik, we obtain the 
eigenstates 

Nk n K + Ns a <?+ 

m= R £ ** - «,)] M |0> « 8 10> * (21) 

whose energy is E — E^ K + En s + E , where E^ K and E^ s are given by 

e n k =puJ2 Vj coth [p(u a - Uj )] - AN K (tt K - N K ) - E K (22) 

ja 

En s = pv £ coth [p(u a - Vj )]- - £ J?- - \ £ J« . (23) 

ja ij i 
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, where E K = J2i + J2j Uij + 2 J2jei s ^i?) . From here on, latin indices, when 
referring to the charge (spin) Hamiltonian, range in Ik (Is), an d greek indices 
from 1 to Nk (Ns). The parameters u> a , vp must fulfill the following set of 
equations 



2 / 

- + X)P cothp[(«j- - u a )] = coth \p(wp - u a )}, (24) 

M 3 

2 / 

- + Z 9 coth ~~ z/ «)] = 2 Z ^ coth [Q( U I3 ~ v o)\ • (25) 

V 3 P 

For XXX- type coupling, present in the spin Hamiltonian (11), that is J z = J, 
the coupling has the form 

■4 ^ iori^j (26) 



Vi - Vj 



the eigenstates are then 

" T ^ - ^ 



^) = nE7T^r.\°)s (27) 



with eigenenergy 



JO! 1 J IJ I 



The quantities z/g must fulfill the RS equations 

- + V^^ = 2V'^^. (29) 



Finally, for uniform Zeeman splitting Q = /i B 5, the eigenstates are still given 
by (27), but the Gaudin equations 

V^^ = 2V'^^. (30) 

j V 3 ~Va p Vp-Va 

have to be fulfilled instead of (29). The eigenenergy E Ns is then 

E "s = E t^zV + vbB(2n s -n s )-W Ji3 (3i) 

ICt " % IJ 1 
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instead of (28), with arbitrary parameters 0j. The couplings are given by 
Eq.(26), replacing Q with 0j. 



4-1 The charge Hamiltonian 



To construct an integrable model for non uniform couplings we generalize 
formula (14) to: 

h k = Y. 2r n f i + E A ij f i f j - E \^ ( 32 ) 

i ij i 



with real symmetric Aij, wherein the operators fj are obtained by substituting 
the isotropic Gaudin models Sj in formula (15) by the anisotropic Gaudin 
Hamiltonians Sj(cr) (see also appendix) 

5>) = £'™>X- (33) 
f^KTf + S,(K). (34) 



A sufficient condition for these operators to commute with each other is that 
the coefficients tug of Ej(cr) can be parametrized as = i;^- = pu coth — Uj)], 



iifp. = yjV. 
l 3 13 



Wij = pu/ sinh [p(ui — tz,-)]. While {wi}, u are real, the parameter 
p can be real or imaginary and switches between hyperbolic and trigonometric 
functions respectively [21]. The resulting Hamiltonian is 



H K = £ 2i*Jf ? - V ^(X+X- + K-K+) + 4 £ (35) 

i ij ij 

+H 3 + H 4 (36) 
with the couplings 

pu{r]i - r}j) 



9ij 



sinh[p(-Uj — Uj)] 



for % 7^ j (37) 
4L/y = + pu(r)i - 773) coth[p(ui - Uj)] , 



The Hamiltonians if 3 and i?4 contain interactions up to tri- and quadri-linear 
in the charge su(2) operators: 
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ff3 = EAj(«?S i (K) + S i (K)A-?) 

ij 

= [A, - A lk ] (v 3k KtK+Kl + w jk K*K*Kl) (38) 

ijk 



^4 = E^A(K)S,(K) 



ijkl ^ 



+ ^ E' - ^ - A^w^K+K- (39) 
+ \ E'E' - ^i* - A n + A kl } WijWki Kf Kj Kf 

ij kl 

+ 7 E'E' {Afe - A,- fc - A« + A kl } VlJ v k iK*K*K z k Kf. 

4 ij fcZ 

It is interesting to note that H% and i/4 vanish for Aij = A, which can be 
also seen directly from equations (38) and (39), because of — (see 

appendix). Another possible simplification can be achieved for A^ = on + oij, 
for which the overall effect of the second term in Eq.(32) is to renormalize the 
parameters r\f 

H K = ^^ + K z a 3 )r y (40) 



We point out that the Coulomb term can be canceled choosing A^ = cosh — uj)] [22]. 
In TCk, the eigenstates of fj(K), and hence of Hk, are 



The eigenvalues ti of f^(K) are given by 

_ ^ # I 

*< = 2 S P COth [P( W a ~~ ~~ 4 Z P COth [P( M i ~~ ~~ 2 " ^ 



The quantities u a must fulfill the equations 
2 / 

- + Ep cothp[(uj - u a )] = 2Ep coth \p(up - u a )] . (43) 
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We discuss the case of constant Ay = A. Then, the Hamiltonian simplifies 
and formula (32) becomes 



H K = Y t 2ru f i + AJ2nf j = Y^2Vif i + A\y2f i \ 

i ij i \ i / 

= Y J ^n + A{K z f. (44) 



In this case Hk is the Hamiltonian (12), provided the coupling constants in 
the latter are parametrized as in Eq. (37) for A^ = A. The eigenstates are 
still given by (41) with the eigenvalues 



E Nk =puJ2 Vj coth \p(u a - Uj)} - AN K {VL K - N K ) - E K . (45) 

jot 

In the p — > limit, which corresponds to the isotropic limit for the Gaudin 
Hamiltonians, the coupling constants become 

9%j = -u(r)i-Vj)/(ui-Uj) , 

for i ^ j (46) 

AUij = A + u(rji - rjj)/ (^ - Uj), 



The eigenstates and eigenvalues are 



i*^=n£-^-io>K 



En k = E Z AN K (Q K - N K ) - E , (47) 

where u) a have to be solutions of the equations 

- + V^^ = 2V'^^. (48) 



Uniform couplings gij = g, Uij = E c are then obtained from a linear u-r] 
relation. Choosing Uj = —2r)j/g, u a = —E a /g, and u — 2, Eqs.(48) are the 
RS equations (17) (we recall that etai = ^ + const.). Setting A = AE C + g as 
well, the Hamiltonian (12) is the BCS Hamiltonian (up to a constant). 
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4-2 The spin Hamiltonian 



We remark, that the procedure leading here to the exact solution of the spin 
Hamiltonian coincides with that applied to H K , because the underlying alge- 
braic structure is still su(2). 
We study the Hamiltonian 



(49) 



with Tj = Sf + Sj(S). Also here, the conservation of the ^-component of the 
total spin restricts the anisotropy of the Gaudin models to the XXZ type. The 
Hamiltonian we get from Eq. (49) is 



ET = £ 2CiSt _ £ ( j ii{S ;sj + sfs]) + J^stsf) , 



(50) 



where the couplings are defined by 

Pv(G - Cj, 



J, 



%3 sinh[p(t>j — Vj)Y 
J*j = -Pv(ti ~ Cj) coth \p(vi - Vj)}, 



(51) 



The eigenstates of H§ xz in Hs and its eigenvalues are: 



iie 



pSt 



sinh \p(is a - vj)] 



l°>S> 



e x n x s =pvJ2 coth [ P (u a - Vj )} - y: q - \ E ^ - \ E 4 

jo. j ij i 

where the u a have to satisfy the RS equations 

-+E — = • 



(52) 
(53) 



(54) 



In order to attack the Hamiltonian (11), i.e. having an XXX- instead of an 
XXZ-type coupling, we have to perform the isotropic limit for the Gaudin 
Hamiltonians and obtain H§ xx = J2i ^d^t ~ Hij -hj Si ■ Sj, with the coupling 
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defined as 



(55) 



The eigenvalues and the eigenstates of Hg xx are 

(56) 

\Y, J h > ( 57 ) 
z i 

For uniform Zeeman splitting, Q = /ie-B (//b is Bohr's magneton, _B the 
magnetic field), the Hamiltonian can be expressed as 

H x s xx = ^bBS z + £(20^.(S) - 3^/4). (58) 

3 

We remark that the constants of motions are Gaudin Hamiltonians themselves 
(instead of the r's). The resulting couplings are given by Eq. (55) with Q = fa. 
The eigenvalues of H s are then 

E Ns ^ P i B B(2Ns-ns)+J2^^-Wj ij -Wj ii , (59) 

ia 1 " Vi 4 ij Z i 

where v a have to be solutions of 

E^^ = 2E'^-— • (60) 

The corresponding eigenstates are given by Eq.(56). The main difference is 
that the v a fulfill the Gaudin Eqs.(60) instead of the RS Eqs.(54). 



5 Discussion of the coupling constants 

In the following we focus on possible choices of the coupling constants in Hk 
that seem appropriate for the physics of metallic grains [23]. It is physically 
reasonable to demand that the pairing attraction decays with the energy dif- 
ference of the interacting levels. This can be achieved by the choice p — 1 in 



^>=nE^-io>5 > 

j a a i j ij 
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(a) 



(b) 



Fig. 1. (a) The pairing and (b) the Coulomb interaction as function of the energy 
difference Ae of the single-particle levels. In Figure (b) the y axis is offset by A/4g 
(see the definition of U in Eqs. (20)) 

Eq.(20) — which corresponds to have hyperbolic functions in the couplings - 
and Ui being a monotonic function of the single particle energies. In order to 
make the couplings ^ and Uij depend on the energy difference Ei—Sj alone, we 
choose Uj = —Ej/Se; 5e has the dimension of an energy (ui are dimensionless), 
and it turns out to be a measure of the range of the pairing interaction. We 
introduce the dimensionless parameter P such that 5e = E^/P, 2Ed being 
the width Ae max of the single particle spectrum. We choose u = Pg/E D) in 
order that the limit P — > corresponds to g^ = g. The resulting pairing is 
plotted in Fig. 1(a). The Coulomb repulsion slowly decreases with the energy 
difference, as shown in Fig. 1(b). For large energy difference and P > 1 the 
decay is linear in the single particle energy difference. 

We now discuss the interpretation of Hamiltonian (8) as a model for M 
weakly coupled grains. We introduce the grain index a = 1, . . . , M (the index 
increasing with the spatial distance of the grains) and the sets I a such that 
for i e I a , Si is an energy level of grain a; U^Li la = I- The charge part e.g. of 
Hamiltonian (8) can then be rewritten as 



-^A/" = X/ X/ £«a C «,iaff Cn . i » ff _ X/ X/ 9i a j b C a,i a ] C a,i a l C b,jbi c b,j b ] 



a=l ia a,b=li a ,jb 
N 

+ E E U iajb n a,iaan b , jb a'- 
a, 6=1 i a ,jb 



(61) 



This Hamiltonian describes a system of many grains, where i a = l,...Q a 
labels the elements of I a and c a ^ aU annihilates an electron with spin a in the 
? a -th level of a-th grain, ^\ For a ^ b, g^ a,b ^ describes the tunneling of Cooper 
pairs; g\ a a ^ = gij if i is the i a -th element of I a , and j the j^-th element of h; 
jj{a,b) j g |- ne Ccmlomb repulsion between the grains a and b. It is written in 
terms of analogously to g^ b \ Couplings g( a > a ) and U^ a ' a ^ describe intra- 
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(a) 



(b) 



Fig. 2. (a) The inter-grain pairing coupling g^ 1,2 \ i.e. the pair-tunneling amplitude, 
between two grains with P = 1 and <E> = $2 — 3>i = 20. The inset shows that the pair 
tunneling amplitude gt 1 ' 2 ) is much smaller than g^ 1 ' 1 ^ = g( 2 ' 2 \ (b) The intra-grain 
Coulomb repulsion is nearly constant for $ 3> P = 1 for fixed pairs of interacting 
levels; this constant increases linearly with Ae (the y axis is offset by A/Ag). 

grain pairing and Coulomb interactions, respectively. Fixing the parameters 
as Uj = $ a — Pej/ E D for j e J a and imposing $ a+i — $ a 3> P, the tunneling 
amplitude of pairs is exponentially suppressed compared to the intra-grain 
pairing coupling, as seen in Fig. 2(a). 



6 Conclusions 



In this paper, we have constructed a new class of integrable models as a func- 
tional of mutually commuting operators, which are a generalization of the 
BCS integrals of motion found in Ref. [12]. The latter are tightly connected 
with isotropic Gaudin Hamiltonians and the generalization presented in this 
paper consists in replacing isotropic with anisotropic Gaudin Hamiltonians. 
The exact solution of the Hamiltonian constructed in this way has been ob- 
tained via the diagonalization of these integrals of motion. This procedure is 
applied to two orthogonal realizations of su(2) (see section 2), leading to the 
model Hamiltonians H K and H s - Eqs. (18) and (19). The resulting models 
describe systems with nonuniform interactions, whose relevance was discussed 
in Ref. [24] according to recent experiments on metallic grains [1]. The nonuni- 
form couplings we obtain are tunable by a set of parameters, which can be 
chosen such that physical demands on the couplings (see Ref. [6,7] and section 
5) are fulfilled. This is demonstrated in Figs. 1 and 2 where the couplings are 
plotted for an exemplary choice of the parameters in the hyperbolic regime, 
p — 1. 

We also discuss possible applications of this class of models to systems of few 
coupled grains and present a choice of parameters such that the intra-grain 
couplings significantly dominate the inter-grain couplings. 
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It is noteworthy that the class of Hamiltonians (18) and (19) can also describe 
a spin chain with a long-range, nonuniform XXZ interaction. Then, r]j and Q 
play the role of a nonuniform magnetic field and the indices (i and j) are to be 
interpreted as labeling lattice sites rather than energy levels. In this case the 
use of trigonometric functions (p = i in Eqs. (20)) can account for boundary 
conditions. 

We want to emphasize that also in nuclear physics [29,30], QCD [31], and as- 
trophysics [32], the BCS Hamiltonian (with uniform couplings) has been thor- 
oughly investigated. The three and four-body interactions, Eq.(38) and (39), 
might come in handy in these contexts. 

Further generalizations of the present work seem interesting. Studies on the 
completely anisotropic XYZ Gaudin model can lead to Hamiltonians where 
the total number of particles (or the total z component of the spin) is not 
conserved. These models are connected with the 2D eight- vertex model [25,26], 
and eventually there is a connection with its off-shell Bethe Ansatz alike that 
between the six- vertex model and the uniform BCS model [14,27]. 
In order to describe interacting bosons, su(l, 1) Gaudin models instead of 
su(2) ones must be considered; for the isotropic su(l, 1) Gaudin model see 
Ref. [28]. Consequently, models based on anisotropic su(l, 1) Gaudin models 
should be studied as well. 

Finally correlation functions for the present class of models can be studied 
along the lines described in Ref. [33] . 



A Gaudin models 

Here we summarize the main results of Ref. [19], which are relevant for this 
work. 

Gaudin studied the following class of Hamiltonians: 



where a a are the Pauli matrices, and the sum over a = x, y, z is implied. 
Mutual commutation of Sj requires 



(A.l) 



k 



(A.2) 



Identifying (x,y,z) with (1,2,3), let us define 



a ■ ea+i( ) ei^.)' 




(A.3) 
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with Uij = —Uji] 1 (m) is an odd function, and Q a+1 (u) are three even func- 
tions. For Uij such that u^ + Ujk + u^i = (e.g. u^ = Ui — Uj), it is found that 
Eqs.(A.2) are satisfied if 6 a are the elliptic Jacobi functions. It is possible to 
find eigenstates and eigenvalues of Sj through Bethe Ansatz [19,20]. 
Imposing conservation of the z-projection of the total spin, X^crf,Sj = 0, 
we get the conditions wj k = w? k , i.e. the anisotropy is of XXZ-type. In this 
case, the (rescaled) elliptic functions reduce to the trigonometric ones 



Wij = w\j = w?j = pw/sinh.[p(ui — Uj)] 



Vij = w^j = pw coth [p{ui — Uj)}. (A. 4) 

Here, we included a real scale factor w together with the parameter p which 
can either be real or imaginary. Defining 

^( u ) = J2 • i r 7 vf (J ? ' 5 " 2 ( M ) = J2P WC ° th iP( U i - M )K> ( A - 5 ) 

i smh[p{Ui -u)\ Y 



the common eigenstates ^Ar T y of Sj for iV| up spins read 

^v T ) = n^ + K*)l°>; ( A -6) 



where |0) is the lowest weight state. The corresponding eigenvalues are 

ii = Y^P W cotn \p( Ui ~~ - 2 J^pw coth [p(tij - uj a )}, (A. 7) 

where u a must satisfy 

cothp[(uj- - u a )] = 2^2' P coth \p(up - u a )]; (A.8) 



a, /? run from 1 to iVj. The isotropic case is obtained imposing [<S, Sj] = 0, and 
corresponds to the limit p — > 0. It leads to the isotropic Gaudin Hamiltonians 



£'^^L (A.9) 



Using here the operators 

^) = E7^L, ^) = Et^i. ( a - 10 ) 
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the common eigenstates of Sj for AT-j- up spins can be written as: 

*iv T )=n^ + K)io), (A.n) 



and the corresponding eigenvalues are: 

^ = E'-^--2E-^- ; (A-12) 



■ Ui - Uj „ Ui-U a 



where u> a have to satisfy 

E^— = 2E'^— (A.13) 

We note, that the specialized Gaudin equations (A.13) and their more gen- 
eral version (A. 8) are the g — > oo limit of the RS equations (17) and their 
generalization (43), respectively. 
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